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1. Introduction 



Meson distribution amplitudes (DAs) describe the momentum fraction distributions 
of partons in a meson, in a particular Fock state, with a fixed number of constituents. 
In the standard treatment of exclusive processes in QCD, which is due to Brodsky 
and Lepage [1], cross sections are expanded in inverse powers of the momentum 
transfer; the size of these power-suppressed corrections, ordered by increasing twist, 
is determined by the convolution of a perturbative hard scattering amplitude with a 
soft nonperturbative DA of given twist. The leading twist DA 0, which describes the 
momentum distribution of the valence quarks in the meson, is related to the meson's 
Bethe-Salpeter wave function (f> BS by 



Here \i denotes the separation scale between perturbative and nonperturbative regime. 
The study of these leading twist 2 DAs has attracted much attention in the litera- 
ture, in particular for the case of the it [2, 3, 4], but only a few investigations are 
devoted to higher twist distributions, which determine the preasymptotic behaviour 
of hard exclusive processes. Higher twist DAs originate from three different sources 
and describe either contributions of "bad" components in the wave function and in 
particular of components with "wrong" spin projection or contributions of transverse 
motion of quarks (antiquarks) in the leading twist components or contributions of 
higher Fock states with additional gluons and/or quark- antiquark pairs. 

DAs of the 7r of twist 3 and 4 have been studied in [5] in the chiral limit, based 
on the techniques of nonlocal operator product expansion and conformal expansion. 
In Refs. [6, 7, 8], vector meson DAs of twist 3 and 4 have been studied, also including 
corrections in the meson-mass. In this paper we extend the analysis of [5] to include 
also terms in the meson-mass in twist 3 and 4 DAs of pseudoscalar octet mesons. 
As discussed in [8], the structure of these mass corrections is more complicated than 
for deep-inelastic lepton-hadron scattering, where the corrections, being induced by 
kinematics, do not involve new information on dynamics and can be absorbed into a 
redefinition of the scaling variable, known as Nachtmann scaling [9]. The situation 
with exclusive decays is different, as matrix elements of operators containing total 
derivatives, specifically 



where is a leading twist operator, vanish for forward-scattering, but do con- 
tribute to exclusive processes. Contributions of the first type can be taken into 
account consistently for all moments of DAs, while contributions of the second type 
are more complicated and can be unravelled only order by order in the conformal 
expansion. Numerically, as expected, these mass terms turn out to be small for the 
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7r, but are dominant for K and the octet DAs of 77. We shall not discuss the octet 
DAs of the 7/ in this paper. The results are of direct relevance for the discussion of, 
for instance, meson transition form factors, 77* — > r/, and also for B meson decays 
into light mesons, see e.g. [10]. 

2. Definition of Distribution Amplitudes 

Amplitudes of light-cone dominated processes involving pseodoscalar mesons can be 
expressed in terms of matrix elements of gauge invariant nonlocal operators sand- 
wiched between the vacuum and the meson state, e.g. a matrix element over a two 
particle operator, 



where T is a generic Dirac matrix structure and we use the notation [x, y] for the 
path-ordered gauge factor along the straight line connecting the points x and y: 



For notational convenience, we refer explicitly to tt~ mesons. For other (nonsinglet) 
mesons, one has to use appropriate SU(3) currents, e.g. l/y/2 (0|-U7^7 5 -u — ^7^75^!^°) 
etc. 

As mentioned in the introduction, we are in particular interested in meson-mass 
corrections. In contrast to vector mesons, whose mass is of order Aqcd and nonvan- 
ishing also in the chiral limit, pseudoscalar meson- masses scale linearly with the sum 
of quark masses. For consistency, we will thus keep all such terms in the analysis 
of the equations of motion, but neglect terms in the difference of quark masses. We 
thus also neglect contributions in the DAs of K mesons that are antisymmetric under 
the exchange of strange and nonstrange quark. 

The asymptotic expansion of exclusive amplitudes in powers of large momen- 
tum transfer corresponds to the expansion of amplitudes like (2.1) in powers of the 
deviation from the light-cone x 2 = 0. As always in quantum field theory, such an 
expansion generates divergences and has to be understood as an operator product 
expansion in terms of renormalized nonlocal operators on the light-cone, whose ma- 
trix elements define meson DAs of increasing twist. To leading logarithmic accuracy, 
the coefficient functions are just taken at tree-level and the distributions have to be 
evaluated at the scale /z 2 ~ x~ 2 . In this section we present the necessary expansions 
and introduce a complete set of meson DAs to twist 4 accuracy. This set is, in fact, 
overcomplete, and different distributions are related to one another via the QCD 
equations of motion, as detailed in Sees. 5 and 6. 




(2.1) 




(2.2) 
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To facilitate the discussion of matrix elements on the light-cone, it is convenient 
to introduce light-like vectors p and z such that 

V, = P, ~ k— , (2-3) 

2 pz 

where P M is the meson momentum, P 2 = m 2 . We also need the projector onto the 
directions orthogonal to p and z: 

gtv = 9ixv-— (p»z„ + p vZll ), (2.4) 
pz 

and will use the notations 

a. = a^, a* = a^jf/ipz), (2.5) 

for an arbitrary Lorentz vector a M . 

We use the standard Bjorken-Drell convention [11] for the metric and the Dirac 
matrices; in particular 75 = i^ ^ 1 ^ 2 ^ 3 , and the Levi-Civita tensor e^ v \ a is defined as 
the totally antisymmetric tensor with e i23 = 1- The covariant derivative is defined as 
D[i = 9^ —igA^. The dual gluon field strength tensor is defined as = ^e^ upa G pa . 

We start with the two particle DAs of the n meson. For the axial vector operator, 
the light-cone expansion to twist 4 accuracy reads: 

(0\u(x) W d(-x)\ir-(P)) = 

= if«P»J o due* Px 4>M+ l -mlx 2 k{u) + % - Um\ ±- x „ J due* Px B(u)(2.6) 

<p n is the leading twist 2 DA, A and B contain contributions from operators of twist 
2, 3 and 4. For brevity, here and below we do not show gauge factors between the 
quark and the antiquark fields; we also use the short-hand notation 

£ = 2«-l. 

The decay constant f w is defined, as usual, as 

(0|«(0)7 M 76d(0)|7r-(P)>=i/ w P /1 . (2.7) 

Numerically, one has f n = 131 MeV and fx = 160 MeV [12]. For i], the situation is 
more complicated due to the mixing with 77', and differing results for the coupling 
to the octet current, are available in the literature. We quote, for instance, 
w 130 MeV [13] and, from a more recent analysis, = 159 MeV [14]. In the 
numerical analysis we will use the ilustrative value = 130 MeV. As we shall see, 
the DAs themselves do not depend critically on the decay constants, although the 
matrix element (2.6) does. 
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The Lorentz invariant amplitude B can be interpreted in terms of meson DAs, 
defined in terms of nonlocal operators at strictly light-like separations, which can 
most conveniently be written using light-cone variables. For the axial vector current, 
the two particle DAs of the ir meson are defined as 

(0\u(z) W d(- Z )\7T-(P)) = 

= if,vJ due^M^ + lUm 2 -^ [ due^ z g w (u). (2.8) 
Jo 1 P z Jo 

Comparing the above with (2.6), one finds 

B(u) = g w (u) - <j> w {u). (2.9) 

The relation of these DAs to those defined by Braun and Filyanov, Ref. [5], is given 
by 

i-tf F {u) = -\ Iim mjB(«), gf F {u)-f dvg* F (v) = -L lim m*A(u). 

Cull Z mi—>0 Jq 10 m£— »0 

(2.10) 

In the local limit — > 0, (2.6) yields the normalization conditions: 



/ 

Jo 



ducp^u) = 1, 

'o 

/ du B(it) = =^ / dug w (u) = 1. 
Jo Jo 

Two more matrix elements define DAs of twist 3 [5]: 

f due* Px <t) p {u), (2.11) 
Jo 



(0\u(x)i l5 d(-x)\ir-{P)) 



m u + m d 

(0\u(x)a af3l5 d(-x)\7r-(P)) = 

i f n ml J f m u + m d x ~ 



. 1 - f ^±^) 1 (P^ - P^) / d«e^^(«).(2.12) 

Also these two DAs are normalized to unity: 

l 

du0(p )ff )(u) = 1. 

The normalization factor in (2.12) differs from the one obtained in [5] by a term of 
0(m u + m d ) ~ 0(ml) : which is tiny for the 7r, but amounts to 10% for the K. This 
may be of particular relevance for calculations of the B — > K decay form factor in 
the framework of QCD sum rules on the light-cone, e.g. [10], to which the twist 3 
DAs give a sizeable contribution. 
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At this point we would like to comment on the numerical values to be used 
for the normalization factors in (2.11) and (2.12). Evidently, it is difficult to give 
precise numbers as long as the quark masses are not more accurately known. To 
circumvent this problem, we invoke chiral perturbation theory (see e.g. [15] for a 
nice introduction), which relates meson to quark masses in the following way: define 
the constant B via the (nonstrange) quark condensate: 

(0\qq\0) = -Y B o 
at the scale /i ~ 1 GeV. Then the meson-masses are given by 

ml = (m u + m d )B , 
m 2 K = (m u4 + m s )B , 

3 1 2 



m 2 m = -\-{m u + m d ) + 2m s \B 0) (2.13) 



where we neglect small corrections in (m d — m u ) 2 . With the standard value of the 
quark condensate, (0|gg|0)(l GeV) = -(0.24 ± 0.01) GeV 3 , one finds B = (1.6 ± 
0.2) GeV. Thus we have 

7T1 2 f 2 



The situation is slightly different for the K and rj (which we consider as a pure octet 
state in this section). Proceeding like with the n, one finds (letting m U (i = 0): 

2 f 2 

m K J K ^ f2 ^ Bo = (Q Q41 ± Q 0Q5) GeV 3 ? 



m s 
m 2 (f 8 ) 2 

m s 3 



(f*) 2 B = (0.045 ± 0.006) GeV 3 , (2.14) 



where one might worry, however, that the constant Bq be affected by SU(3) violation. 
Using the actual values for the meson masses and m s (l GeV) = 150 MeV, one finds 

m 2 f 2 m 2 (f 8 ) 2 

= 0.042 GeV 3 , W = 0.042 GeV 3 , 
m s m s 

which is in good agreement with the results from chiral perturbation theory. 

Let us now define the three particle DAs. To twist 3 accuracy, there is only one: 

{0\u{z)a^gG aP {vz)d{-z)\'K'{P)) = 
f m 2 

= i — (p a p^ - p a p u g^ - V0Pv.9v* + VnVvgi^ T{v,pz) + ... (2.15) 

77T>u ~r TTld 
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where the ellipses stand for Lorentz structures of twist 5 and higher and where we 
used the following short-hand notation for the integral defining the three particle 
DA: 

T(v,pz) = J ' Vae- ipz ^- ad+va ^T(a d ,a u ,a g ). (2.16) 

Here a is the set of three momentum fractions ad (d quark), a u (u quark) and a g 
(gluon). The integration measure is defined as 

J Da — J dadda u da g 5(l — a u — ad — a g ). 

There are also four three particle DAs of twist 4, defined as 

(0\u(z)w 5 gG af3 (vz)d(-z)\n-(P)) = 

= PuiPaZp - VpZa) — f n mlA\\(v,pz) + ippg^ - Pa9p ll )fnmlA±(y, pz)(2.17) 
pz 

(0\u(z)^igG af3 (vz)d(-z)\7r-(P)) = 

= P^(PaZ(3 ~ Pf3Z a ) — f w mlV\\(v,pz) + (ppgap ~ Pa^)/ W "l^V±(u, pz) . (2.18) 

pz 

A short synopsis of the various light-cone projections of the three-particle matrix 
elements and their relation to DAs is given in Table I. 



Twist 


(/ii/aft) 










3 


■ ± ■ ± 


T 
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■ ■ * 


A\\ 


V|| 








_L± • 


A ± 





Table 1: Identification of three-particle DAs with projections onto different light-cone 
components of nonlocal operators. For example, _L_L • refers to V>7_l75Gj_.?/'- 



For completeness, let us mention that also four particle twist 4 DAs exist, cor- 
responding to contributions of Fock states with two gluons or an additional qq pair. 
Such distributions will not be considered in this paper for two reasons: first, it is 
well known [16] that four particle twist 4 operators do not allow the factorization 
of vacuum condensates such as (ipip), (G 2 ). Because of this, their matrix elements 
cannot be estimated reliably by existing methods (e.g. QCD sum rules), although 
they are generally expected to be small. Second, and more importantly, the four 
particle distributions decouple from the QCD equations of motion in the two lowest 
conformal partial waves. To this accuracy, therefore, it is consistent to put them to 
zero. Vice versa, nonvanishing four particle distributions necessitate the inclusion 
of higher conformal spin corrections to distributions with less particles, which are 
beyond the approximation adopted in this paper. 
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3. Conformal Partial Wave Expansion and Equations of Mo- 
tion 



The aim of this paper is to express the DAs defined in the previous section in a 
model-independent way by a minimal number of nonperturbative parameters. The 
one key ingredient in solving this task is the use of the QCD equations of motion 
which will allow us to reveal interrelations between the different DAs of a given 
twist. Nonlocal operators on or near the light-cone can conveniently be treated in 
the framework of the string-operator technique developped by Balitskii and Braun 
[17]. In the present context, we need the following nonlocal operator identities [8]: 

d 



dx 



- u(x)-f^ 5 d(-x) = -i / dvvu(x)x a gG aiJ ,(yx)'y fl 'y 5 d(-x) 
ti J-i 

+ (m u - m d )u(x)ij 5 d(-x), (3.1) 
d p {u{x)^ p ^d(-x)} = -i J dvu(x)x a gG aix (vx)-f^ 5 d(-x) 

+ (m u + m d )u(x)i^ 5 d(-x), (3.2) 

d f 1 

d, l u{x)a fJ , v ^d(-x) = -i 7^- u(x)j 5 d(-x) + J dvvu(x)x p gG pv (vx)j 5 d(-x) 

-i J dvu(x)x p gG plx (vx)a IJiU ^ 5 d(-x) 
+ {m d - m u )M(x)7 i ,7 5 rf(-a:), (3.3) 
- u(x)a flv j 5 d(—x) — —id v u(x)j 5 d(—x) + / dvu(x)x p gG pu (vx)^^d(-x) 

IX J-i 

-ij dvvu(x)x p gG ptl (vx)a flv j 5 d(-x) 
- (m u + m d )u(x)^ 5 d(-x). (3.4) 



d_ 

dx 



Here d p is the total derivative defined as 



d 

d p {u(x)Td(-x)} = — {u(x + y)[x + y, -x + y]Td(-x + y)} 



By taking matrix elements of the above relations between the vacuum und the 7r~ 
meson state, one obtains exact integral representations for those DAs that are not 
dynamically independent. 

The other key ingredient in our approach is the use of conformal expansion [18, 5] 
which, analogously to partial wave decomposition in quantum mechanics, allows one 
to separate transverse and longitudinal variables in the wave function. The depen- 
dence on transverse coordinates is represented as scale-dependence of the relevant 
operators and is governed by renormalization-group equations, the dependence on 
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the longitudinal momentum fraction is described in terms of irreducible representa- 
tions of the corresponding symmetry group, the collinear conformal group SL(2,R). 
The conformal partial wave expansion is explicitly consistent with the equations of 
motion since the latter are not renormalized. The expansion thus makes maximum 
use of the symmetry of the theory in order to simplify the dynamics. 

To construct the conformal expansion for an arbitrary multi-particle distribution, 
one first has to decompose each constituent field into components with fixed Lorentz 
spin projection onto the light-cone. Each such component has conformal spin 

3 = 2 + 

where / is the canonical dimension and s the (Lorentz) spin projection. In partic- 
ular, / = 3/2 for quarks and / = 2 for gluons. The quark field is decomposed as 
ifj + = (l/2)^ip and ip_ = (1/2)^/0 with spin projections s = +1/2 and s = —1/2, 
respectively. For the gluon field strength there are three possibilities: G.± corre- 
sponds to s = +1, to s = —1 and both Gj_± and G.* correspond to s = 0. 

Multi-particle states built of fields with definite Lorentz spin projection can be 
expanded in irreducible unitary representations of SL(2,R) with increasing conformal 
spin. The explicit expression for the DA of a m-particle state with the lowest possible 
conformal spin j — ji + ■ ■ ■ + j m , the so-called asymptotic DA, is 

rh (rv rv rv \ - T< ^ ]l j L 2jm ^ rv 2 -?' 1 " 1 rv 2 -?' 2 " 1 rv 2 ^ 1 

PasKOi!, a 2 , , a m ) — — — a l a 2 ■■■a m ■ [o.Q) 

T(2j 1 ) ■ ■■r(2j m ) 

Here are the corresponding momentum fractions. This state is nondegenerate and 
cannot mix with other states because of conformal symmetry. Multi-particle irre- 
ducible representations with higher spin j +n,n = 1, 2, . . ., are given by polynomials 
of m variables (with the constraint Y^k=i a k — 1 )> which are orthogonal over the 
weight-function (3.5). 

In particular, for the leading twist 2 DA 0^ defined in (2.6), the expansion goes 
in Gegenbauer polynomials: 

0.(m,/i 2 ) =6u(l-u) ^l + fja 2n (/i 2 )C 2 3 / t 2 (2 M -l)j . (3.6) 

To leading logarithmic accuracy, the (nonperturbative) Gegenbauer moments a n 
renormalize multiplicatively with 

a n {Q 2 )=L~<~ /b a n {^) ) 

where L = a s (Q 2 )/a s (fj, 2 ), b = (lliV c — 2Nf)/3, and the anomalous dimension 7„ is 
given by 

^ = 4Cf (" ( " + 2)+ ^-I- 2( n+ l)(n + 2) 
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In this paper, we work to next-to-leading order in conformal spin and thus truncate 1 
the above expansion of (f) n after the term in n — 1. For the ir, the corresponding 
Gegenbauer moment was determined in e.g. [3] from QCD sum rules, for K, we use 
the value determined in [10]: 

oT 2 {\ GeV) = 0.44, of (1 GeV) = 0.2, aP 2 {\ GeV) = 0.2. (3.7) 

The value for rj is new and follows from an analysis of the QCD sum rule in [3] by 
fixing the continuum threshold sq to reproduce /® = 130 MeV. 



4. Meson-Mass Corrections 



The structure of meson-mass corrections in inclusive processes is in general more 
complicated than that of target-mass corrections in deep inelastic scattering, which 
can be resummed using the Nachtmann variable [9]. The terms entering the Nacht- 
mann variable are just the subtracted traces of the leading twist forward-scattering 
matrix element, which can also for exclusive processes be summed to all orders. Let 
us illustrate this point with the aid of the two-point correlation function of scalar 
fields: 

1 



(0\c/>(x)<f>(-x)\P) 



due iiPx 



ip(u) + - x 2 m 2 ip 2 (u) + 0(x 4 ) 



Here ip is the leading twist DA; ip2 receives contributions from both the subtraction 
of traces in the leading twist operator, also dubbed "kinematical" corrections, and 
from intrinsic "dynamical" higher twist corrections. The subtraction of traces can 
be done on the operator level by making use of the condition [17] 

d 2 



[cf>(x)<f>(-x)] lx = 0, 



(4.1) 



dx Q dec q 

which translates into the condition that all local operators arising in the Taylor 
expansion be traceless. A formal solution is [17] 

[cf)(x)(f)(-x)} lt = 

d 2 



= 4>(x)(f)(-x) + 



CO 

E 

n=l 



dtt 



1 , 
— x 1 
4 



[t(l-t)\ 



n-1 



t Qx Q doc q 



4>(tx)4>(—tx). 



To order x 2 , one thus has 

(f)(x)(j)(-x) = 

= [<l>(x)<i>(-x)] u 



1 f 1 

-x 2 / dttd 2 [(j)(tx)(f)(— tx)]\ t + interaction terms +0(x 4 ). 
4 J o ' (4.2) 



1 Note that a thorough discussion of the shape of the w DA of leading twist necessitates the 
inclusion of higher terms in the conformal expansion. In this paper, however, we concentrate on 
higher twist DAs which constitute corrections to the leading twist DAs that are suppressed by 
powers of the characteristic momentum transfer in hard reactions, so we feel justified in neglecting 
higher order conformal corrections to these corrections. 



9 



Note that d 2 [. . .]i. t . is a higher twist operator. Taking a forward-scattering matrix 
element, (P\ . . . |P), the second term on the right-hand side vanishes, and all mass 
corrections arise from subtracting traces in the leading twist matrix element. This 
is the source of the Nachtmann corrections. In the following we will calculate the 
corresponding matrix element for exclusive processes, 

(O|[#r)0(-z)] u .|P), (4.3) 

exactly, i.e. summing up all terms in x 2 . 

It is clear that (4.3) can only depend on the leading twist DA ip(u). Writing 
(4.3) as a Taylor-series over the moments of ip, we can subtract the traces explicitly 
for each term in the expansion, which yields (P 2 = m 2 ): 

(O|[0(x)0(-x)] u .|P> =X) i j «°»» U mV U n (Px/V^). (4.4) 

n=0 n ' V / 

Here U n are the Chebyshev polynomials 



j=0 

and ((O n )) are the moments of ip: 



[f] . _ A 

o— n \ J J 



((On)) = f 
JO 



duCi>{u). 



Using the generating function of U n , it proves possible to sum (4.4) explicitly (see 
also [19]): 

(O|[0(x)0(-x)] Lt .|P> = 

= / duV>(u)— == — e W2 sin(f a/C^P - m 2 x 2 /2) . (4.5) 

Jo y/(Px) 2 - m 2 x 2 du L vv ; / ;j v ; 

In the spirit of Nachtmann, we would like to absorb all terms in m 2 into a new scaling 
variable Px — > (Px)'. Although we did not succeed in finding such a variable, the 
above expression can be simplified considerably by introducing 

2{Px)' = Px+ ^(Px) 2 - m 2 x 2 , 

so that 

(o\[<Kx)<K-x)] Lt .\p) = 

= f 1 du ( eX p{^(Px)'} - exp li^P x y^L^-\) .(4.6) 
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Expanding to 0(x 2 ), this can be written as 



(O|[0(x)0(-x)] Lt .|P) = J dt I (1„ 
and combining with (4.2), we get 



e* Px + -m 2 a?t£ 2 e* Pxt 
4 



(O|0(x)0(-a;)]|P) 



dt / du 



This means that both sources of mass corrections, the subtraction of traces in the 
leading twist matrix element and the higher twist operator containing total deriva- 
tives, act in the same direction and thus enlarge the mass correction terms. We will 
observe the same effect, enlarged mass corrections, also in QCD. 



5. Twist 3 Distribution Amplitudes 

The twist 3 DAs of the tc have already been studied in [5]. Here we extend this study 
by including terms in p\ = (m u + rrid) 2 /m 2 , ~ 0(m 2 ). 

To next-to-leading order in conformal spin, the only three particle DA T gets 
expanded as 

T(a) = 360r] 3 a u a d a 2 g j 1 + u 3 ^ (7a g - 3) j . (5.1) 

rj 3 is defined as 

(tt\ua^gG a pd\-K~) = 

2 

Tfl 

= ifnV3 ~ ~ P a Pu9^ ~ PfiP^va + PpPuQa^) , (5.2) 



and ui 3 is defined as 



3 

(Q\wT IJt Z'y 5 [iD li , gG a ^]d - - idpua^gG^d^') = 

f m 2 3 

= % — — 2P a P^P II — V3 u 3 + 0(higher twist). (5.3) 

m u + ma Zo 

In the notations of [5]: 

These parameters are scale-dependent with (Ca = N c ) 

m (Q 2 ) = L^%(n% l v 3 = fc F + C A , 

25 7 

lu 3 (Q 2 ) = L^usOm 2 ), 7s =-6^+3 C * 
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Numerical values are obtained from QCD sum rules [5, 20] and collected in Table 2. 

The two particle DAs PjCr are determined by T and p\4>^. As an analysis of the 
matrix elements of the exact operator relations in Sec. 3 leads to integro-differential 
equations that cannot be solved in a closed form, we prefer to perform an analysis 
of moments: 

Ml = 5 n0 + ^ Ml_ 2 + 2(n - l)Mll 2 + 2(n = 2) M?_ 3 - M*_ 2 , 

Ml = 5 n0 + M°_ 2 + m* + -pl-^- Ml (5.4) 

n+3 n 2 n+3 ni n+3 nl w n+3 n v 7 

Here we use the notation M£ = J* du£ n (f)p(u) and the functions 

Tl = / c?«d / da u — T(a), 0r 2 = / da d / da u — (a d - a u - ^) T(a). 
Jo Jo « 3 Jo Jo "a 

Except for the new terms in p^, the relations for moments agree with those obtained 
in [5]. 

Conformal expansion imposes that 4> p gets expanded in Gegenbauer polynomials 
ClJ 2 and <j) a in Cn" 2 [5]. From the recursion relations for moments we find: 

M«) = i + (so^ - \ pI) c^io + (-3^-3 pI^) cl /2 (0, 

<t> a {u) = 6u(l - «) {l + (5% - \ ifeo* - ^ Pi - \ Pla^j C 3 2 /2 (0 1 . (5.5) 



In Fig. 1 we plot the two two particle DAs of twist 3 for ir, K and 77 mesons. 
Evidently, the effect of mass corrections is not negligible and for <p p even modifies the 
shape of the DA near the endpoints, which is due to the dependence of the coefficient 
of cT on p\. We would like to recall, however, that the above parametrizations are 




Figure 1: The two particle DAs of twist 3: (f> p (u) (a) and 4> a {u) (b). 
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to be understood in the sense of (mathematical) distributions rather than as models 
that are valid point by point, and that they are always intended to be convoluted 
with smooth perturbative scattering amplitudes, which in particular will smooth out 
the effect of neglected higher order terms in the conformal expansion. 





7T K 1] 


a 2 

V3 

fi 2 [GeV] 2 


0.44 0.2 0.2 
0.015 0.015 0.013 

-3 -3 -3 
0.0077 0.096 0.12 



6. Twist 4 Distribution Ampli- 
tudes 

In this section we repeat the analysis of 
twist 4 DAs performed in [5] in a more sys- 
tematic way and extend it by including mass 
correction terms. Table 2: Input parameters for twist 3 

Due to G-parity, in the chiral limit, the DAs > calculated from QCD sum rules. 
DAs .A,, and A± are antisymmetric under The accuracy is about 30%. Renormal- 
the exchange of a d and a u , whereas Vy and ization scale is 1 GeV " 
Vj_ are symmetric; contributions of "wrong" G-parity give rise to asymmetric con- 
tributions to the two particle DAs of K and are neglected in the following. The 
distributions A\\ and V\\ correspond to the light-cone projection 7.G.* (see Table 1) 
and have the conformal expansion 

V\\(a) = 120a u a d a g (v 00 + v 10 (3a g - 1) + . . .), 

A\\(a) = 120a u a d a g (0 + a 10 {a d - a u ) + . . .), (6.1) 

respectively. Note that the leading spin contribution to A\\ vanishes because of G- 
parity (for massless quarks). 

The DAs V± and A±, on the other hand, correspond to the projection 7j_G'._i_ 
and thus do not describe states with a definite projection of the quark spins onto 
the light-ray z^. We separate the different quark spin projections with the aid of the 
auxiliary amplitudes T~0^ and defined as 

(0\u(z)igG a /3(vz)'y.'y li 'y*d(-z)\ir-) = f n ml (ppg^ - p a g^) H n (v,pz), 

(0\u{z)igG a p{vz)^^.d{-z)\n-) = j^m\ (ppg^ - p a g^) H l \v,pz). (6.2) 

The original distributions A±_ and Vj_ are then given by 



V ± («) = --(H n («)+^ T («)), 

A ± (a)= l -{H n (a)-H l \a)). 
and have a simple expansion in terms of Appell polynomials, to wit: 



(6.3) 



rOHa) 



n u \a) = 60a d a 2 , 



60a u a g 



hoo + h 01 (a g - 3a d ) + h w ( a g - 



hoo + h i(a g — 3a u ) 



-a, 



hw [ a g - -a d 



(6.4) 
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where we have taken into account the symmetry properties, i.e. 

H n {a d , a u ) = H l ^(a u , a d ). 
From (6.3), the following relations can be derived immediately: 



V±(a) = -30at 



h 00 (l- a g ) + V 



a g (l - a g ) - Qa u a d 



+h 
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<Xg(l-<*g) ~ 2^ + a d) 



A±(a) = 30a 2 (a u - a d ) 



h o + h 01 a g + - h w (5a g - 3) 



(6.5) 



The DAs V± y \\ and A±,\\ depend, to next-to- leading accuracy in the conformal spin, 
on a total of six parameters: t>oo an d ^oo of leading conformal spin and t>io, aio, 
h 10 and h Q1 of NLO conformal spin. Our next task is to relate these parameters to 
independent local matrix elements. Defining 



{0\u^igG^ a d\n ) = f n mlr] 4 P a , 



(6.6) 



which is equivalent to 



(OlwyJgG^dl-K ) = -- Uml^P^g^ - P v g m ), 



it follows 



^oo = ^oo = -g Va- 



(6.7) 



774 is scale-dependent with 



r? 4 (Q 2 )=L^/W), H = \C F . 

To NLO in conformal spin, beyond the matrix elements already defined above, 
we need only one more matrix element of a conformal quark-gluon operator: 



(0\u[iD^igG v ^d- -id^uigG^d^ ) = 
= f^ml^uJi [p^P u - J mlg^ + O (twist 5). 
The scale-dependence of w 4 is given by 

u A (Q 2 )=L^/ b uj A (fi 2 ), 74 



(6.8) 



10 r 8 



Numerical values for 774 and 0J4 were calculated from QCD sum rules [21, 5] and are 
collected in Table 3. In the notation of Ref. [5], 5 2 = m 2 7]4, e = 21/8 cu 4 . 
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The procedure how to relate i>io, aio, h w and h i to local matrix elements is 
described in detail in Ref. [8], so that we mention only the essentials. Three of the 
four necessary relations follow from an analysis of various light-cone projections of 
the matrix elements of the operators 



° = Dp gG pu + gG pv i Dphad, 
°2U = 5/3 gG^ + gG pv i D^cH^ 

The fourth relation can be derived from the operator identity 



\pn E mf$ = -12iu7 P 75 \g pP D a - D a G p/3 + (a <-> /?)} d - Ad p u(-i^G ap + -f a Gp p )d 

8 ~ 8 ~ 28 ~ 

- g dpu-y a G aa d - - d a u~i a G afi d + y g a pd p u-i^ 5 G ap d, 



where 

"15 - - , 3 



«7m75 £> Qj D/3 d - -dadpwy^d - traces 



- symmetrized 

is a leading twist 2 conformal operator. The matrix element of this operator is pro- 
portional to the Gegenbauer moment a2 of the twist 2 DA and brings in a dynamical 
mass correction in the twist 4 DAs. This is precisely the effect mentioned in the 
introduction: total derivatives of twist 2 operators enter higher twist DAs. 
The final results for the NLO parameters read: 

21 9 
21 

VVQ = — 7?4^4, 
o 

7 3 

V = 4^4 - 7^02, 

7 3 

hio = 2 ^ 4a;4 + 20 ° 2 - ( 6 - 9 ) 

For a 2 — >■ 0, these results agree with those obtained in [5]. 

We are now in the position to derive expressions also for the remaining two DAs 
of twist 4, g n and A. From the operator relations (3.1) and (3.2), one obtains 

d f u f u 2 
g n (u) = 2(p p (u) - <p n (u) + — / da d da u — (A\\(a) - 2A±(a)) , (6.10) 

« M Jo Jo a g 

A(«) = 12/ dv dw{g v (w)-<f> v (w))-2 dv(2v - l)(<f>M + g*(v)$.U) 
Jo Jo Jo 

+ / da d / da u — (a d - a u - £)(2A±(a) - A\\(a)). (6.12) 

Jo Jo a a 
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g w corresponds to a definite quark spin projection and thus has a simple expansion 

1/2 

in Gegenbauer polynomials C n ■ 

oo 

9n (u)^J2^ C 2i 2 (0- (6-13) 

i=0 

From (6.10), one finds: 

9o = 1, 

18 20 

92 = 1 + y ^2 + 60% + y 774, 

where we neglect terms of 0{m u + ma) induced by <f> p , as g n itself enters the matrix 
element of the axialvector current already as 0(m 2 ) = 0(m u + ma), cf. (2.8). 
The expansion of A is not that straightforward and involves logarithms: 

. . _ f 16 24 20 

A(m) = 6mm < — + — a 2 + 20% + — % 

+ {~h + Te ~ h ™ ~ i *) + ("^ - " 4 

+ (-y «2 + 21%^ {2m 3 (10 - 15m + 6m 2 ) hiM + 2m 3 (10 - 15m + 6m 2 ) InM 

+mm(2 + 13mm)} . 



In Fig. 2 we plot m 2 p gp and mpA for the mesons P = ir, K, rj. Whereas the DA 
mpgp is not too different for the three different mesons, the impact of meson-mass 




Figure 2: The two particle DAs of twist 4: gp(u) (a) and A (b) for P = ir, K, g. 



16 



corrections on mpA is more noteworthy. For the area under the curve, i.e. the overall 
normalization Np = dum P A(u), we find: 

N n = 0.47, N K = 0.70, N v = 0.77. 

The result for N v is essentially independent of the precise value of /®. The impact 
of meson-mass corrections is thus rather profond for the DA A. Likewise the change 
in normalization of <p a , this shift in A may have a noticeable effect on the B — > K 
form factors calculated from QCD sum rules on the light-cone. 

7. Summary and Conclusions 

In this paper we have studied the twist 3 and 4 two 
and three particle DAs of pseudoscalar nonsinglet 
mesons in QCD and expressed them in a model- 
independent way by a minimal number of nonper- 
turbative parameters. The work presented here is 
an extension of the paper [5] on 7r DAs and is in 
particular concerned with corrections in the meson 
mass. The one ingredient in our approach is the use 
of the QCD equations of motion, which allow one 
to express dynamically dependent DAs in terms of 
independent ones. The other ingredient is conformal expansion which makes it pos- 
sible to separate transverse and longitudinal variables in the wave functions, the 
former ones being governed by renormalization group equations, the latter ones be- 
ing described in terms of irreducible representations of the corresponding symmetry 
group. 

The analysis of twist 4 DAs is complicated by the fact that the twist 4 terms are 
of different origin: there are, first, "intrinsic" twist 4 corrections from matrix elements 
of twist 4 operators. There are, second, admixtures of matrix elements of twist 3 
operators, as the counting of twist in terms of "good" and "bad" projections on 
light-cone coordinates does not exactly match the definition of twist as "dimension 
minus spin" of an operator. There are, third, meson-mass corrections, which one 
may term kinematical corrections, that come, on the one hand, from the subtraction 
of traces in the leading twist operators and, on the other hand, from higher twist 
operators containing total derivatives of twist 2 operators. Meson-mass corrections of 
the first kind are formally analogous to Nachtmann corrections in inclusive processes, 
while the contribution of operators with total derivatives is a specific new feature 
in exclusive processes, which makes the structure of these corrections much more 
complex. 

Our final results are collected in Sees. 5 and 6. In contrast to the case of vector 
mesons whose mass is finite also in the chiral limit, pseudoscalar octet mesons bring 





7T 


K 


V 


a 2 


0.44 


0.2 


0.2 


V3 


0.015 


0.015 


0.013 




-3 


-3 


-3 


Va 


10 


0.6 


0.5 


UJ4 


0.2 


0.2 


0.2 



Table 3: Input parameters for 
twist 4 DAs, calculated from 
QCD sum rules. The accuracy 
is about 30%. Renormalization 
scale is 1 GeV. 



17 



in the complication that their mass squared depends linearly on the quark masses, so 
that for consistency one also has to take into account terms of 0((m u + m^) 2 /m%) ~ 
0(ml). The effect of meson-mass corrections is noticeable for twist 3 DAs and in 
particular for the twist 4 DA A, whose normalization increases by ~ 50%, when 
comparing the n with the K. 

We hope that our results will contribute to a better understanding of SU(3) 
breaking effects in hard exclusive processes and in particular to the investigation of 
B and B s decay form factors into n, K and r\ mesons from QCD sum rules on the 
light-cone. 
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